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We consider the inequivalent quantizations of a iV-body rational Calogero model with a Coulomb 
(3JT). type interaction. It is shown that for certain range of the coupling constants, this system admits 

a one-parameter family of self-adjoint extensions. We analyze both the bound and scattering state 
| sectors and find novel solutions of this model. We also find the ladder operators for this system, 

with which the previously known solutions can be constructed. 
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D ' 1. INTRODUCTION 

Exactly solvable quantum many body systems like Calogero model and its variants [H 0, 0| have found diverse 
applications in many branches of contemporary physics, including generalized exclusion statistics [4], quantum hall 
effect Tomonaga-Luttinger liquid quantum chaos 0], quantum electric transport in mesoscopic system 
^2 1 spin-chain models , Seiberg-Witten theory [11] and black holes [12j . The rational Calogero model is described 

^vq . by N identical particles interacting with each other through a long-range inverse-square and harmonic interaction 
ly-^ [ on the line 1]. The exact spectrum of this rational Calogero model with harmonic confinement has been found 
through a variety of different techniques [l|, [3| , all of which impose the boundary condition that the wavefunction 



OO 
O 



and the current vanish when any two or more particles coincide. With this boundary condition the Hamiltonian 
is self-adjoint, which ensures the reality of eigenvalues as well as the completeness of the states. However it was 
. £^ I found later that, within a certain region of the parameter space, there exist more general boundary conditions for 
which the rational Calogero Hamiltonian (with and without harmonic confinement) admits self-adjoint extensions 



and yields a rich variety of spectra [171 118| | . As is well known, the possible boundary conditions for an operator are 



encoded in the choice of its domains, which are classified by the self-adjoint extensions [l9| of the operator. Such 
self-adjoint extensions play important roles in a variety of physical contexts including Aharonov-Bohm effect 2C|, 
two and three dimensional delta function potentials 2l| . anyons 22| , anomalies 23( , (^-function renormalization 2J] , 



particle statistics in one dimension [25| and black holes [26j . So it should be interesting to find out more examples of 



exactly solvable models which can be quantized by using the method of self-adjoint extension. 

In this context it may be noted that an exactly solvable variant of the rational Calogero model has been constructed 



by Khare 27j, where the confining simple harmonic potential is replaced by a coulomb- like interaction. The bound 



states of this model can be related to those of the rational Calogero model with harmonic confinement by using the 
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underlying SU(X, 1) algebra [28[. However, apart from having an infinite number of bound states, this model with 
coulomb-like interaction also supports continuous scattering states. Similar to the case of original Calogero model, 
these bound as well as scattering states have been constructed by using the boundary condition that wavefunction 
and the current vanish when any two or more particles coincide. Due to a factorization property of the eigenfunctions, 
the eigenvalue problem of this many body system can be reduced to that of the corresponding radial Hamiltonian. 
In this article, our aim is to find out more general boundary conditions which admit self-adjoint extensions for the 
radial part of the rational Calogero model with coulomb-like interaction and study the related spectra. 

The arrangement of this article is as follows. In Sec. 2 we briefly recapitulate how to separate the radial part of 
the rational Calogero Hamiltonian with coulomb-like interaction. Then we discuss about the most general form of 
eigenstates associated with this radial Hamiltonian H r . These eigenstates could be singular or nonsingular at r = 



value of the radial coordinate r. The bound and scattering states found by Khare are all nonsingular at r = [27|. In 
Sec. 3 we show that such nonsingular bound states can also be constructed by using creation annihilation operators 
associated with the underlying SU(1, 1) algebra. In Sec. 4 we show that the radial Hamiltonian H r admits self-adjoint 
extensions within a certain region of the parameter space. Inequivalent quantizations of H r by using this method lead 
to bound and scattering eigenstates which are in general singular at r = 0. We explicitly construct such bound states 
and scattering states in Sec. 5 and Sec. 6 respectively. In Sec. 6, we also derive the scattering matrix for the scattering 
states and show that the eigenvalues of the bound states can be reproduced from the poles of this scattering matrix. 
Sec. 7 is the concluding section. 

2. GENERAL FORM OF EIGENFUNCTIONS OF THE RADIAL HAMILTONIAN 

The Hamiltonian for the rational Calogero model with the Coulomb type term is given by 

N q2 1 



, dx: 

1=1 L i<Cj 



It describes the TV-body problem with equal mass in 1-dimension and in units such that 2m/i~ 2 = 1. It is also 
understood that the coupling constant g is constrained to satisfy condition g > — i. The parameter a will be 
allowed to have any value, positive or negative, which will provide us with the possibility to treat the attractive as well 
as the re puls ive Coulomb potential, respectively. This point is somewhat different in comparison to the treatment 
made in 2JJ where only the case for a > is considered. Having the Hamiltonian ([I}, we intend to solve the 
eigenvalue problem 

= (2) 

Following [l[ , we consider the above eigenvalue equation in a sector of configuration space corresponding to a definite 
ordering of particles given by x\ > X2 > ■ ■ ■ > Xn- The translationaly invariant eigenfunctions of the Hamiltonian 
H can be factorized as 



* = l[(x i -x j ) a+ ^<f>(r)P k (x), (3) 



\X% Xj 
i<j 

where x is an abbreviation for {x%, X2, • • •, xjy}, a — ±i^/l + 2g and r is the collective radial variable defined as 

= -^£(si-3*) a . (4) 



r 2 

N ■ 

i<3 



Functions Pk(x) are translationaly invariant, homogeneous polynomials of degree k, k > 0. They satisfy the 
equation 



E^ + ^+te^ 



d 2 1 v-^ 1 d d 

^ 4 Pk(x) = 0. (5) 



dxf 2 (xi — Xj) dxi dxj 

i=l 1 ?:<j J 
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This equation is analyzed in detail by Calogero [lj]. After inserting the factorized form of from Eq. Q into Eq.([2|) 
and using Eqs.Q and (O, we get an equation for <j)(r) as 

H r <f>= E(j), (6) 

where 

ft + 1,1!-^. (7) 

The parameter b which enters Eq.([7]) is defined as 

, N(N-1), 1 N AT 3 

It should be noted that, due to Eq. {3]), the measure of the quadratically integrable <j)(r) is given by da = r°dr, 
where f3 = 2k + 2b + 1. Thus, the eigenvalue equation ([2|) of the Calogero model with Coulomb- like interaction is 
reduced to the eigenvalue equation (J6j> of the corresponding radial Hamiltonian H r , which may be explicitly written 
as 

- -^0(r) - (2k + 2b+ l)~4>(r) - -?=-<Kr) = Ecb(r). (9) 
ar z r ar J Nr 



However, for the purpose of finding out all possible boundary conditions on the wavefunctions for which H r admits 
self-adjoint extensions, in due course we shall also need to study the deficiency subspaces for this radial Hamiltonian. 
Thus for the sake of convenience we consider a somewhat more general form of Eq. ^ , namely, 

- -^(r) - (2k + 2b+ l)~<Kr) - -^=-0(r) = E^r), (10) 
dr 2 - r dr JNr 



with E — E,+i,—i, depending on whether we are interested in the eigenfunctions or the deficiency subspaces for 
the radial Hamiltonian ([7]). 

By performing the following transformations 

(f>(r) =r~%ip(r), /3 = 2fc + 26 + 1, y = cr, (11) 



with the parameter c as yet unspecified constant, Eq. (|10p can be reduced to Whittaker's equation, which in turn is 
related to confluent hypergeometric equation of the form 

y^ + (2M + i-y)^-(5 + /.-«)]x(») = o. (12) 

The last step, where Whittaker's equation, satisfied by the function tfi, reduces to Eq. (TT2")) , is performed by means 
of the factorization ip(r) = e~%yix(y)- The parameters introduced in Eq. (fT2|) are given as 

f3 1 

c = 2\^E, (13) 



cVN y/-4NE 

The general solution to Eq. p2D is known to be the linear combination of the confluent hypergeometric functions of 
the first and the second kind [29j 

X (r) = AM (1+h-k, 2/x+l, cr) +BU C~+ p- K, 2/i+l, cr) , (14) 



4 



where A and B are arbitrary constants. Due to the relations ip(r) = e ^y^x(y) an d V = cr : Eg. (p~4|) yields the 
general form of ip(r) as 

i;(r) = e-^(cr)^ ^AM ^ - n, 0, crj +BU - k, /?, cr^j . (15) 
Since the measure on the space of quadratically integrable <j> functions is da = r^dr, Eq. Jlip implies that we have 

4>*4)d(j = / ip*tj)dr, 



showing that on the space of quadratically integrable ip functions the measure is simply dr. 

It may be noted that, in terms of some generic parameters a, b and the variable z, the confluent hypergeometric 
functions M and U are determined by the expressions [2^ 

az a(a + l)z 2 (a)„z™ 

M ^^^ 1 + T + ^Trk + --- + T^! + ---' (16) 



7T 

U (a, b, z) 



M(a,b,z) 1 _ b M{\ + a-b,2~b,z) 



Y{l + a-b)Y{b) v ' T(a)T(2-b) 



(17) 



sin irb 

where the symbol (a)„ means 

(a) n = a{a + l)(a + 2) (a + n - 1), (a) = 1. (18) 

These expressions show that while M(a, b, z) is nonsingular at z = 0, U(a, b, z) could be singular at z — with leading 
power z 1_b . Thus, for B ^ 0, the solution of ip(r) given in Eq. (|15p is singular at r = 0. In this article, our main aim is 
to employ such singular solutions of ip(r) to construct bound and scattering states through the method of self-adjoint 
extension. On the other hand it may be observed that, the solution of ip(r) given in Eq. (|15l) would be nonsingular at 
r = for the case B = and (3 > 0. Bound and scattering state solutions found in Ref. [27J all correspond to such 
nonsingular solutions of ^(r). 

3. CONSTRUCTION OF NONSINGULAR BOUND STATES THROUGH LADDER OPERATORS 

Let us now concentrate on the set of bound states corresponding to the spectrum found by Khare. These states 
are solutions to Eq. lfTO]) when E = E, E < 0. In this case, parameters c and n, introduced in (fT3)) . become real 



parameters c = 2^/|i?| and n = —j==. Since U is singular at r = 0, after utilizing equation (fTTj) . we are left with 
the wavefunctions of the form 

<j>{r) = r-h~^ cr (cr)^M(^ - k, (3, cr). (19) 

Due to the fact that M comprises an infinite diverging series, it has to be truncated and this can be achieved by 
setting 2 ~ K = ~ n ' n = 0)1)2,3,... With this truncation condition M reduces to the associated Laguerre 
polynomials, L„ (a;) = ^ , " M (— n, a + 1, x), and wavefunctions (fH))) belong to the following discrete set of 
normalized eigenfunctions, labeled by the quantum number n, 



Mr) = y (2n + /3 " r(ri + /3) c^e-^L(f- 1 )(cr) = C^-^L^ier), n = 0,1,2,3, ... (20) 

They are normalized to unity with respect to the measure da = r^dr, i.e. they satisfy J (f>*<f> n da = 1. The 
corresponding bound state energies follow from the aforementioned truncation condition and are given as 

E n = ~^r - f 175> n = 0,l,2,3,... (21) 

4A/ (k + b + n + i) 
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By using recursive relations for the associated Laguerre functions 

(n + a)4-i W = (2n + 1 + a - ,)LW(x) - (n+ 



(22) 



we can find recursions 



x^L£\x) + (n + l + a- x)L^(x) = (n + ljL^x), 



1 d\ n(2n + /3-2)(n + /3-l) ^ 
" - 2 V - V Ty) MV) = V 2^3 



(23) 



! (2n + P + 2)(n + P)(n+ 1) 



2n + P 

satisfied by the radial functions l|2"0|) . If we introduce the number operator defined as 

Ncf> n = n^> nj 



<?Wi(y), 



we can easily find ladder operators 15(] , 30(] , 3l|] from the recursive relations (|23|) and (pM|) . They are 

b = 



(2/V + /3) 



(iV+/3-l)(2JV + /3-2) 



(24) 



(25) 



(26) 



6t 



(27V + /?) 



(/V + /3)(2/V + /? + 2) 

Straightforward calculation shows that the ladder operators 6 and are bosonic, 

[6, & + ] = 1, 

together with 

[iV, &] = -&, [JV",6 + ] = 6 + , 



(27) 



(28) 



(29) 



resulting in the simple relation including the number operator TV = b^b. If we take the vacuum state to be the 
wavefunction (I20p with the lowest possible energy (which happens when n = 0), it is seen that b annihilates this 
vacuum state, namely, b(f>o = 0. The bosonic operators b and b^ satisfy the usual oscillator relations, 



i>n = Vn4> n -x, b'4> n = V n + 1 </> n+ i. 



(30) 



Due to relations (|10p . (f2Tj) and (|25|) . the radial Hamiltonian in Eq. ([7]) can be expressed in terms of the number 
operator N in the way 



H r = — 



4N {N+k+b+^y 

We would like to factorize the radial Hamiltonian (|3ip. so that it can be written in the form 

Hr = -LlU + fc, 

7 



(31) 



(32) 
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with 7 and h as yet undetermined parameters. In order to make such factorization, we have to make a transition 
from the pair of bosonic oscillators {b, tf} to new pair of operators {A, A^}. Although the form of the Hamiltonian 
(1321) looks much more simple, the price we have payed is that the new deformed oscillators are no more bosonic, but 
rather they are deformed, obeying a more complicated relation in terms of the number operator, namely, 



A^A = ®(N). 



(33) 



In the above relation $(iV) is an analytic function which is required to satisfy following three conditions: 

(«) <P(N) > 0, (34) 



(m) $(0) = 0, 



(35) 



(iii) $(1) = 1. 



(36) 



The function $(iV) which is consistent with relations (|32[) and (f33|) and which obeys conditions (i),(ii) and (iii) is 
given as 



$(A>) 



7 2 a 2 / 1 



1 



AN \ (P 



(f) (tf+f) 



(37) 



The parameters h and 7 are introduced so as to accommodate for the condition (ii) and the normalization condition 
(iii), respectively, and are equal to 



Nf3 2 



(3\% + 1) N 
a 2 (f3+l) ' 



(38) 



The deformed oscillators A, A^ can be related 32| to the bosonic oscillators (|26p and (f2"T)) in the following way 



V N V N 

If we further introduce the operators J+ , J_ , Jo defined as 



J_ = A 



N 



(39) 



(40) 



'$(7V) \J®{N) 

one can show that they are, in fact, generators of SU(1, 1) algebra, 

[J_,J+] = 2J , [J ,J±]=±J±. (41) 

In papers [3], the underlying conformal symmetry of the rational Calogero model with a Coulomb-like term 
is revealed by constructing an explicit realizations of the corresponding SU(1, 1) generators. These realizations 
happen to be different from those found in 1J], l_5j where the realizations of SU(1, 1) generators for the rational 



Calogero model with the harmonic confining term are considered. Since it is known 14|, 15[ that all models with 
underlying conformal symmetry can be mapped to the set of decoupled oscillators, one could do the same for the 
Hamiltonian |T]) by using the construction of SU(1, 1) generators made in [28J], [33|. After finding an appropriate 



similarity transformation, one could apply it to SU(1, 1) generators to find ladder operators for the Hamiltonian 
It is possible to carry out such transformation since the all systems with underlying conformal symmetry have radial 
excitations described by the associated Laguerre polynomials with two of the generators playing the role of creation 
and annihilation operators in the equivalent problem including decoupled set of oscillators. This approach would lead 
to ladder operators which would not coincide with the ladder operators (|40"]) . but would rather be related to them by 
means of some similarity transformation. 
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4. DEFICIENCY INDICES OF THE RADIAL HAMILTONIAN 

The spectrum of this model discussed above is valid for the usual boundary conditions where the wave function 
vanishes at r = and it is square integrable. We shall now find the most general set of boundary conditions for which 
the radial Hamiltonian H r is self- adjoint. For this we follow the method of von Neumann. We start by recalling the 
essential features of this method 



19] 



Let T be an unbounded differential operator acting on a Hilbert space TL and let D(T) be the domain of T. The 
inner product of two element a, (3 G TL is denoted by (a,/3). Let D(T*) be the set of <fi E TL for which there is a 
unique r) e TL with (T£, <j>) = (£,77) V £ G D(T). For each such 4> G D(T*), we define T*(f> = r\. T* then defines the 
adjoint of the operator T and D{T*) is the corresponding domain of the adjoint. The operator T is called symmetric 
or Hcrmitian iff [T(j>, rj) = (<p, Trj) V (j), T) G D(T). The operator T is called self-adjoint iff T = T* and D(T) = D(T*). 

We now state the criterion to determine if a symmetric operator T is self-adjoint. For this purpose let us define 
the deficiency subspaces K± = Ker(i =p T*) and the deficiency indices n±(T) = dim[i£±]. Then T falls in one of the 
following categories: 

1) T is (essentially) self-adjoint iff (n+,n_) = (0,0). 

2) T has self-adjoint extensions iff n + = There is a one-to-one correspondence between self-adjoint extensions of 
T and unitary maps from K + into K_ . 

3) If n + =^ ri_, then T has no self-adjoint extensions. 

We now return to the discussion of the effective Hamiltonian H r . This is an unbounded differential operator defined 
in R + . H r is a symmetric operator on the domain 

D(H r ) ee {0(0) = <f>'(0) = 0, 0, (/)' absolutely continuous, </> G L 2 (cZct)}, 

where da = r@dr. We would next like to determine if H r is self-adjoint in the domain D(H r ). To perform such an 
analysis it is necessary to obtain the square-integrable solutions of the equation 

H*cj> ± {r) =±#±(r). (42) 

The operator H* is the adjoint of H r and is given by the same differential operator as H r , although their domains 
might be different. Below we shall give the analysis for the parameter range where \x ^ 0, the case for /1 = being 
similar. Thus, Eq. (|42| is identical to Eq. (fTU|) when E = ±i, 

- -T2<t>±{r) - (2* + 2b+l)~cf> ± (r) - -^=-£t(r) = ±#±(r). (43) 



We are interested in finding the square-integrable solutions to Eq. (|4"5|) . The solutions of Eq. (|4"2")) or Eq. (T53J) which are 
square-integrable at infinity are given by <f>±(r) — r~T?p±(r), where 

V>±(r) - e-^± r (c ± r)^C/(^ - k±, /?, c±r) (44) 

with 



(45) 




Since these solutions are also required to be square-integrable near the origin, it is necessary to investigate their 
behaviour for r — » 0, which looks as 1 



g 7T 

ip±[r) — > (c±r) " 



sin 7r/3 



1 (c^) 1 ^ 



r(i - 1 - «±)ro9) r(|-«±)r(2-/3) 



(46) 



1 In subsequent considerations we shall work on the space of ip(r) functions where the measure is dr. 
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In the above expression we have restricted ourselves to the lowest order in r, so that we could take M(a, b, z) — ► 1 
as the argument z tends to 0. The square integrability of the wavefunction (|44[) near the origin is determined by 
(HHJ) which implies that as r — » 0, 



\ip±(r)\ dr 



A ir f3 + A 2 r + A 3 r 2 ~P 



dr, (47) 



where A\, A%, A% are some constants independent of r. From Ea.(|47p it is seen that near the origin, the functions 
ip± (and consequently functions <j>± ) are not square- integrable for the parameter (3 satisfying f3 < — 1 or (3 > 3. 
Consequently, in the parameter range (3 < — 1 or (3 > 3, the functions ip± are not the elements of the vector space 
L 2 [R + ,dr] of quadratically integrable functions defined on the positive real axis. In that case, n+ = n_ = and 
H r is essentially self-adjoint in the domain T>(H r ). However, if — 1 < (3 < 3, the functions ip± (and consequently 
functions <j>± ) are square-integrable. Thus, if (3 lies in this range, we have n + = n_ = 1 and the Hamiltonian H r 
is not self-adjoint in the domain T>(H r ), but admits self-adjoint extensions. Note that from (|13p. the allowed range 
of (3 implies that the parameter jj, must lie in the range —1 < fi < 1. 

The above allowed range of [i, together with |5]), ([TT]) and (fTB")) , implies that the values of N, k and a + | must 
satisfy the relation 

iV-l + 2fc 1 iV-5 + 2fc 

- N(N-l) <a+ 2 < - N(N-l) - (48) 

for the self-adjoint extension to exist. For N > 3, we have the following classifications of the boundary conditions 
depending on the value of the parameter a + \ . 



(i) a + | > | : This corresponds to the boundary condition considered by Khare in [27J. For this choice, both the 
wave-function and the current vanish In this case, /j, > 1 for all values of k > 0. The corresponding 
Hamiltonian is essentially self-adjoint in the domain D(H r ), leading to a unique quantum theory. 

(ii) < a + i < i : For this choice we see that the wave-function in ([3]) vanishes in the limit Xi Xj , although 
the current may be divergent. In this case fj, > and k must be equal to zero so that n may belong to the range 
< fx < 1. The corresponding constraint on a + 5 is given by < a + | < jyfj^n ; which can only be satisfied for 
2V = 3 and 4. So new quantum states associated with the self-adjoint extension of H r exist only in the k = sector 
of N = 3 and N = 4. 

(iii) — i<a+i<0: The lower bound on a + h is obtained from the condition that the wavefunction be square- 
integrable. The parameter a + ^ in this range leads to a singularity in the wavefunction 'J in Eq. ([3|) resulting from 
the coincidence of any two or more particles. Using permutation symmetry, such an eigenfunction can be extended 
to the whole of configuration space, although not in a smooth fashion. The new quantum states in this case exist for 
arbitrary N and even for non-zero values of k. In fact, imposing the condition that the upper bound on a + i should 
be greater than —i, we find from (|48)l that k is restricted as k < \ (iV 2 — 3./V + 10). It can also be shown that there 
are only two allowed values of k when both N and a + 5 are kept fixed. 

Von Neumann's method also provides a prescription for obtaining the domain of self-adjointness of a symmetric 
operator, which admits a self-adjoint extension. The extended domain T> z (H r ) in which H r is self-adjoint contains 
all the elements of T>(H r ), together with the elements of the form e l iip+ + e~ l ^ip~, where z € R (mod 2ir). Thus 
the self-adjoint extensions of this model exist when — 1 < f3 < 3, and in that case, 

V z (H r ) = V(H r ) © {e l ^+ + e -*fy_} 

is the extended domain in which H r is self- adjoint. 



5. BOUND STATES OF THE RADIAL HAMILTONIAN WITH SELF-ADJOINT EXTENSION 



We shall now find solutions of the physical problem for the range of system parameters where the self-adjoint 
extension is necessary. In finding the solutions to Eq.(frD|), we shall first consider the bound state sector of the problem. 
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In this sector, the energy of the system is negative, E < 0, and the wavefunctions need to be square-integrable. We 
consider the solution of the form 



ip(r) = Be~^(cr) 2 U 



(3 



k, /3, cr 



(49) 



In order to make an analysis and to find the spectrum, we should know the behaviour of the U function near the 
origin. Using Eqs. (fT6| and (Tl~7|l . we can expand U(a,b,z) at z — > limit as 



U(a, b, z) 



sin nb 



1 



f 



r(i + a-&)r(&) 
i 



a a(a + 1) z 2 
1+ b Z + b(b+l)2! 



0(z 3 



2-b 



r(a)r(2 - b) 

Consequently, at r — > limit, ip(r) in Eq. (|49|) behaves as 

1 



a _ ^ + (l + a-fr)(2 + a-6) z 2 + 3 



(2-6)(3-&) 



2! 



(50) 



V>(r) — ► Be-5 cr (cr) 



£ 7T 



sin 7r/3 



(cr) 



r(i-f 

1-/3 



«)r(/3) 



l 



cr 



r(f - «)r(2 - p) 



i 



i 



2-/3 



K (1 

— cr H 



(f -«)(§-« + !) 
2/903+1) 

-4_«)(2-|-«) 



2(2-/3) (3-/3) 



(cr) 2 + 0(r 3 ) 



(cr) 2 + 0(r 3 ) 



(51) 



The parameters c and k appearing in (|5ip are given in ||13|) . except that they are evaluated for E = E. Since the 
energy is negative, these parameters are real, 

c = 2-J-E = 2\/K= Pl 



(52) 



CVN y/-4NE y/WE b py/N 

Here, for future convenience, we have introduced the absolute value Ef, of the bound state energy, Ef, — —E, E < 0, 
and the real parameter p which coincides with c in the bound state sector. Note that c and k will no more be 
real in the scattering sector. 

If the wavefunction (|49p is expected to describe a physically acceptable bound state solutions to Eq. ([1^) , it has to 
belong to the domain of self-adjointness T> z (H r ). If ipo(r) £ D(H r ), then an arbitrary element of the domain T> z (H r ) 
can be written as ipo( r ) + p{^ l ^^P+ + e ~ 1 ^ ip-)t where p is a constant. If the solution of the physical wavefunction |49|) 
belongs to the domain T> z (H r ) , the functional form of physical wavefunction must match with that of an arbitrary 
element of the domain T> z (H r ) , which is given by 



-0W = ipo{r) + p(e l *ip + + e 



(53) 



Inserting Eqs. (|5ip and (|46|) into relation (|53|). and equating the coefficients of the lowest order powers in r (for which 
there is no contribution from ip (r)), yields the following two conditions 



Bel 



e. 

e 2 c+ 2 



2 C_ 2 



r(i-K-f) r(i -«+-§) r(i-K. 



2 / 



e 2c+ 



r(f-K) r(f-« + ) r(fi-«_)' 

where B = -B/p. After dividing both sides of these two expressions, we get the relation 



(54) 



r(i-«-f) „i- g _ r ( #-. + ) 



r(f-K-) 



r(f-K) 



/3 



(55) 



r(i-«+-"fj + r(i-K_-#) 
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Inserting the expressions (|52|) for c and k we obtain the final condition 

@ ot e 1 ^ C+ 1 , e~ z ^ 



2 s/SWEj' (r, fWY p r(f-« + ) r(f-«_) 
— — ^ - (2yg b ) - | — | , (56) 

r(i-« + -fj r(i-«_-f) 

which determines the spectrum corresponding to bound states of the radial Hamiltonian H r {7} as well as the initial 
many-body Hamiltonian (jTJ) . 

Writing — -0 = £ie t01 and c+ — 3- = &e 102 , ([5*6]) can be expressed as 



r(|-^_) (2 ^ } = U + (57) 

The above analysis shows that for a given choice of the system parameters, Eq. ([57]) gives the energy eigenvalue 
E = —Eb as a function of the self-adjoint extension parameter z. For a fixed set of system parameters, different 
choices of z lead to inequivalent quantization and to the spectrum for this model in the parameter range where the 
system admits self-adjoint extension. In general, the energy E — —Ef, cannot be calculated analytically and has to be 
obtained numerically by plotting ([57]) . Figures 1 and 2 show l.h.s and r.h.s of Eq. ([57|) for two different, representative 
sets of the system parameters as well as for the two different choices of the self-adjoint extension parameter z. The 
curved lines at those figures represent graph of the function f(Eb) which is given by the l.h.s of Eg. (|57|) . On the other 
hand r.h.s of Eq. (|57p is represented by a horizontal straight line. The energy eigenvalues of the system described by 
the Hamiltonian (JTJ) are obtained by looking at the intersections of these two curves. We see from figures that there 
is an infinite number of bound states near Eb — > 0. For a > 0, there are infinite number of bound states for any 
value of z. However, the existence of non-oscillatory part shows that, just like the usual case, the spectrum has a 
lower bound for all possible values of z. The situation when a > is shown at figures 1 and 2. 

For the choice of the self-adjoint extension parameter z = z\ such that 61 + 4^ = f , the r.h.s. of (f57j) is zero. This 
implies that 

which gives the usual energy eigenvalues as expressed in (f2Tj) . It can be shown that the choice of z = z-i such that 
O2 + 4f = § gives a similar result. At this point it may be noted that the analytical solution ([58]) implies that for a 
certain values of the self-adjoint extension parameter and system parameters, even the repulsive Coulomb potential 
leads to the formation of only one bound state. It can easily be seen if we write (f551) in the form a = \/NEb(2n + (3). 
This expression shows that in order to have the repulsive Coulomb potential, that is a < 0, one has to restrict (3 
within the range — 1 < (3 < and set n equal to zero, resulting in a single bound state. The same conclusion holds 
also in the general case where the analytical solution is not possible, and it can be verified by extensive numerical 
investigation of the general relation (|56[) (see Figure 3 as an example). 

Finally, it is important to emphasize that the system described by |T]) has a fundamentally different behaviour 
depending on the sign of a. While for a > 0, the l.h.s. of Eq. (f5T]) exhibits oscillatory as well as non-oscillatory 
behaviour, leading to infinite number of bound states, for a < 0, it shows only non-oscillatory behaviour resulting 
in the existence of at most one bound state. This single bound state, if it exists, shows up only for the certain range 
of the self-adjoint extension parameter z. For a = 0, this observation is consistent with the result obtained in [16| . 
This feature can most easily be seen by looking at the special case (58) where the analytical solution is available. 

There, in order for a to be less than zero, we must have n — together with (3 within the range — 1 < (3 < 0, 

2 

resulting in a single bound state E = —Eb = —jfm, as already stated just after Eq.(58). 
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Figure 1. A plot of Eq. (f57|) using Mathematica with N = 1000, a = 50, (3 = 0.8, k = 1 and z = 0.1. The horizontal straight 
line corresponds to the value of the r.h.s of Eq, (|57[) . 



f (E b ) 




Figure 2. A plot of Eq. ([57]) using Mathematica with N = 100, a = 1.5, j3 = -0.7, k = 1 and z = -0.73. The horizontal 
straight line corresponds to the value of the r.h.s of Eq, (|57|l . 




Figure 3. A plot of Eq. (|57[) using Mathematica with iV = 1000, a = —1, j3 = 1.5, k = 1 and z = 0.1. The horizontal straight 
line corresponds to the value of the r.h.s of Eq. (|57p . This graph shows the general feature exhibited for arbitrarily strong 
repulsive Coulomb potential, i.e. for any a < 0. 
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SCATTERING STATES OF THE RADIAL HAMILTONIAN WITH SELF-ADJOINT EXTENSION 



Let us now turn our attention to the scattering sector of the problem described by Eq.JTU]). Since the scattering 
states correspond to positive energy solutions of Eq. (flQ|) when E = E > 0, the variable y — cr = 2\/—Er becomes 
purely imaginary, i.e. y = iqr, where the real parameter q is defined as q — 2y~E . Therefore, for analyzing the 
r — ► oo limit of the scattering states, it is of importance to know the behaviour of the confluent hypergeometric 
functions M(a,b,z) and U(a 1 b,z) in the asymptotic region Re(z) = and Im(z) — > +oo. Following Abramowitz & 
Stegun, we can expand confluent hypergeometric functions in this asymptotic region as 



M(a,b,z)~^^\e z z a - b 
T(a) 



1 + 0(1; 



T(b - a) 



(-*)" 



l + 0{\z\- 1 ) 



(59) 



U(a,b,z)^0(\z\- a ). 



(60) 



Due to the fact that we are dealing with the problem where Re(z) = 0, both leading terms in the asymptotic 
expansion (|59[) of M approximately have the contribution of the same order, so that both of them have to be taken 
into account. In order to find the scattering matrix we could equally well take the following linear combination 



X {y) = AM Q + M - K , 1 + 2 M , y^j + By~^M Q - fx - k, 1 - 2/i, y^j 



(61) 



as a general solution to Eq. (TT2"|) . instead of the one given in ([14")) . In this case the solution for the function tjj, appearing 
in (TTTj) would look like 



r(r! = , 3" (cr )l ( A (q)M (| - «, f3, cr) + B{q){crf- p M (l - | - k, 2-/3, 



(62) 



where we have assumed that the coefficients A(q) and B{q) depend on the real parameter q. 
By using (|59p . we have the following r — > oo limits 



M(| - k, P, cr) 



r(|-«) 



e cr (cr) 



T(f3) , 



r(f + «) 



(63) 



M(l - | - k, 2 - /3, cr) 



r(i-f-«) 



r(2-/3) 
r(i-f + «) 



(-cr) 



(64) 



so that the wave function 



describing the scattering state, in the above limit behaves as 



ip(r) = ip(E = E) — ► A(q) 



r(f-K) 



"(cry K +A(q) 



r(f + «) 



e 2' 



r (2-/3) „i 



r(i-f -«) 



e* or (cr) _ ' , + B(«j)(-l) 



r(i-f + «) 



(cr) 



(65) 



Note that the parameter k is also purely imaginary in the scattering sector. For the coupling constant a greater 
than zero, k can be expressed as k — —i ^-^ = —i j^L = — i\n\. By using the relations y = cr = iqr and k = — 
we can express ip( r ) i n Eq. (|65p in terms of oscillatory incoming wave and outgoing wave as 



ip(r) = ip(E = E) 



A(q) 



r(4-«) 



s(?) 



r(2 - p) 
r(i-f -«) 



„i(^5T+|it| lnr) i 



e l z K q K 



A (q)^- e M K -Z) 

v w r(§ + «) 



ff(V) F ( 2 e i7r( K +|-l)| e -i(i«r+|K|lnr) 

r(i-| + K ) i 



(66) 
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The scattering matrix and the corresponding phase shift can be obtained from the above limiting form of the wave 
function as a ratio of its outgoing and incoming amplitudes, 

S(a) = e v( * q) = 2 — 2 (67) 

Next, to find a relationship between so far unspecified constants A(q) and B(q), we use the expansion (I16p to obtain 
the r — > limit of the wave function (1621) . in the lowest order in r, 



tp(E = E) > A{q){cr)~- + B{q){cr) 1 ~^. (68) 

We recall that the Hamiltonian H r admits a self-adjoint extension in the parameter range 3 > /3 > — 1. Since the wave 
function (|6"2"|) has to belong to the domain of self-adjointness T> z (H r ) — T>(H r ) © {e l '-0 + + e -1 ^ V-} we can write 

pi}j{E = E) = e'^V+ + e~ if V-, (69) 

where p is some constant and, as before, ip± are square integrable solutions of Eq. (|10|) when E = ±i, respectively. In 
the limit r — > 0, the behaviour of ip± is given by the relation (|46|) . Since according to Eq . (|69|) . the coefficients of 
appropriate powers of r in (|68p and (|46|) must match, the following two conditions emerge 

p4Mc I= e 5 - ^ A" 1 he '5- - R , 70) 

smTT^r^- § -K+)T(fi) sin7r/3r(l-|-/t_)r(/3) 

pS^c 1 ^-^ — ! e-^— C — (71) 

sin7r£r(f -K+)T(2-p) sin^/3 r(| - «_)r(2 - /?) ' 

The last two equations yield 

i. §. 
Ml) T(2 - P) e ' f r(i-f-* + ) + e ~'* r(i-gltt-) ^ 

5(g) r(/3) iA e+1 -# , _ ia c _i-4 c • 1 ' 

r(f-/s+; r(f-K_) 
By using this expression, the scattering matrix (|6T[) becomes 

5(g) = e 2 *^) = FlUi ' a ' Z) r(l 7 } g' 2K , (73) 

*i(£> «, z) r(f+K) r(l-f +k) 

where c = 1\J —E = ig, and ft = = ^ 2ne ' ^ n wri ting the expression for the scattering matrix we have 

introduced the following two functions 

F ' (Aa ' 2> = e ' f it^ + e "' 5 i^Z)' (74) 

a(ft a. ,) = e-i r(i _Y_ ^ + .-t r(i ^ ■ (75) 

where z is the self-adjoint extension parameter and c± and /c± are defined in (|45|) . As a remark, one can note 
that the functions F\ and i*2 are simply related as 



F 2 ((3, a, z)=F 1 (2-0, a, z). 
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As it is seen from the form of the scattering matrix, for any given value of (3 in the parameter range which admits 
self-adjoint extension, the scattering matrix has an infinite set of poles on the positive imaginary axis of the complex 
g-plane. The existence of poles for the scattering matrix means that there are bound states in the system under 
consideration. By taking q = ip as some arbitrary pole for the scattering matrix (173|) . one can obtain the following 

2 

equation determining the bound state energies Ef, = —E = ^- : 

F 2 (J3, a, z) e ^^-P+h) q l -P e «r(#+«-X) ^ 

Fi(/3, a, z) r(| + k) r(i-f + k) ~ ' ( } 

This expression, after utilizing the set of relations q = 2\f~E = ip = i2\fE~l and k = — i-^^ = = ~ 2 sfW V/v ' 

finally gives 

~ 2 ~ 2^VW^ 1-/3 Fi(A a, z) 
which reproduces the bound state condition (I55|) . 



CONCLUSIONS 



In this paper we have analyzed the TV-body rational Calogero model with a Coulomb like interaction. We have 
shown that for certain ranges of the system parameters, the system admits a one parameter family of self-adjoint 
extensions. The results obtained here for both bound and scattering state sectors are very different from those 
obtained by Khare in 27). However, there is no contradiction between these findings as they refer to different ranges 
of the system parameters. We have also shown that for specific choices of the self-adjoint extension parameter, the 
usual results of Khare can be recovered. 

It has also been shown that a ladder operator construction exists for this system, which also leads to the solution 
found by Khare. This construction indicates that su(l,l) can be regarded as a spectrum generating algebra for this 
system, as it happens in conformal quantum mechanics 3J] yielding equispaced energy levels. We think that there 
is a strong correlation between our and the constructions made in papers 28|. 331 ] - We hope to address this issue in 
more detail in a future. 

In the presence of the self-adjoint extension, the su(l,l) can no longer be implemented as the spectrum generating 
algebra as the dilatation generator in this case does not in general leave the domain of the Hamiltonian invariant 
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21 



23l [35|. As a result, the spectrum for a generic choice of the self-adjoint extension parameter is no longer 
expressed in the Coulomb-like form. However, when z — z\ or z%, the Coulomb- like nature of the spectrum is recovered 
and su(l,l) can again be implemented as a spectrum genera ting algebra. This effect is analogous to the quantum 
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3G] 



anomaly also observed in the pure Calogero type systems 

We have also seen that the system exhibits qualitatively different behaviour on two sides of the point a = 0. We 
find that for the attractive Coulomb potential (a > 0) there exists an infinite number of bound states. In the case 
of the repulsive Coulomb potential (a < 0), there appears to be at most a single bound state, which exists only for 
certain values of the self-adjoint extension parameter. 

In this paper we have restricted our discussion to the case when the coupling constant g of the inverse square 
interaction is such that there is no collapse to the centre. It would be interesting to analyze this problem where the 



coupling is more attractive with g < — |, which would require renormalization group techniques 2ll. |37| 
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